Slow light in three-level cold atoms: a numerical analysis 
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We investigate theoretically the slow group velocity of a 
pulse probe laser propagating through a cold sample and in- 
teracting with atoms in a three-level A configuration having 
losses towards external states. The EIT phenomenon pro- 
duces very small group velocities for the probe pulse in pres- 
ence of a strong coupling field even in presence of the pop- 
ulation losses, as in an open three-level system. The group 
velocity and the transmission of the pulses are examined nu- 
merically as functions of several parameters, the adiabatic 
transfer, the loss rate, the modification of the atomic velocity 
produced within the cold sample. The conditions for a more 
efficient pulse transmission through the cold atomic sample 
are specified. 

PACS: 42.50.Gy, 32.80.-t, 42.50.Vk 



I. INTRODUCTION 

In a a three level A system, with a central excited 
level connected by electric dipole transitions to the two 
ground ones, electromagnetically-induced transparency 
(EIT) and reduction in group velocity are based on 
the low-frequency coherence created by the laser radi- 
ation [0|2j. One optical transition is driven by a strong 
laser, denoted as coupling laser, and the second tran- 
sition is driven by a weak laser, the probe one. The 
probe laser absorption and dispersion are determined by 
the modifications in the populations and in the coher- 
ences produced by the coupling laser, and also by the 
low-frequency coherence between the ground states cre- 
ated by the simultaneous application of the coupling and 
probe lasers. For the three-level atomic A configuration, 
the absorption is decreased by the coherent trapping of 
the population in the ground state [§. The probe fight 
group velocity is greatly reduced by the very steep dis- 
persion associated with the coherent population trapping 
narrow resonance. Remarkable results for the reduction 
of the group velocity in three-level systems were recently 
reported in cold samples as well as in samples with 
thermal velocity distributions 

Often the real atomic transitions of quantum optics 
phenomena do not correspond to the simple theoretical 
models, and several multilevel atomic structures can be 
described as open systems, i.e., systems where popula- 
tion, whence coherent population trapping, is lost be- 
cause of decay into sink levels not excited by the lasers. 



Even if the open features of a multilevel system may be 
eliminated through the application of a repumping radia- 
tion to deplete the sink state , they play an important 
role in most experiments. In this work we plan to inves- 
tigate the propagation of slow light in an open atomic 
system for the specific conditions of the cold atom ex- 
periment by Lau et al jij. We analyze a sample of cold 
sodium atoms laser excited in an open A scheme starting 
from two different hyperfine levels and with losses from 
the excited state towards another ground state. The ex- 
periment of 0] was performed with both a cold sodium 
atomic sample and a sodium condensate. However the 
atomic interactions are not included in our analysis [pT| . 

The present work addresses several issues related to 
the slow light propagation. Our main aim is to deter- 
mine the role of the open three-level system losses on the 
group velocity, probe pulse amplitude transmission and 
transmission bandwidth. In an open three-level system, 
any excited state occupation leads to a loss of the atomic 
population towards atomic states not excited by the laser 
light, whence to a modification of the slow light propaga- 
tion. Our numerical simulations determine the laser pa- 
rameters more appropriate to realize both efficient pulse 
slowing and transmission. We have specifically investi- 
gated if the STIRAP configuration could be used 
to reduce the role of the excited state occupation in an 
open three-level system. In STIRAP the counterintuitive 
coupling/probe pulse sequence with detuned lasers pro- 
duces a very efficient coherent-state atomic preparation 
from the initial ground state to the final ground state, 
and decreases the role of the excited state occupation 
with losses towards external levels. We have verified that, 
even if the counterintuitive pulse sequence enhances the 
probe pulse transmission, the EIT conditions for slow 
light do not correspond precisely to those for the STI- 
RAP process. We investigate also the influence of the 
adiabatic transfer conditions on the pulse transmis- 
sion: a larger probe pulse amplitude transmission is ob- 
tained for a larger temporal width of the probe pulse, 
because that width modifies the adiabaticity condition. 

The transmission of a slow probe pulse depends on the 
decay rate of the coherence between the lower levels of 
the A system. An efficient pulse propagation, i.e. a small 
absorption coefficient, is obtained by reducing the coher- 
ence decay rate. In a sample of room temperature atoms, 
collisions and transit time are the main contributions to 
that decay rate. We examine the influence of a coher- 
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ence relaxation rate on the pulse transmission and on 
the velocity reduction. In a cold atom sample, the previ- 
ous contributions to the coherence decay rate are greatly 
decreased. However in a cold sample with long interac- 
tion times the kinetic energy motion of the atoms rep- 
resents another important contribution to the coherence 
decay rate. The forces by the laser on the atoms modify 
the atomic momentum and whence the parameters of the 
laser atom interaction. We investigate the pulse propaga- 
tion taking into account the atomic momentum modifica- 
tions occurring in the laser interaction. The importance 
of light forces on atoms in the case of slow light propaga- 
tion has been recently investigated jl5| for the configu- 
ration where a slowly propagating laser pulse produces a 
force acting on a second atomic species contained within 
the propagation medium. In the present work the light 
forces act on the same atoms producing the slow light 
process. Thus we investigate the limitations imposed on 
the slow light propagation by the light forces acting on 
the momentum of the atoms producing the light slowing. 

The present analysis makes use of the time/space de- 
pendent numerical solution the Optical Bloch Equations 
(OBE) of an open three-level system. Section II intro- 
duces the three-level system and the coupling/probe laser 
propagation. Section III presents the numerical results 
for the propagation of a laser pulse through a cold sample 
of open three-level atoms in the presence of a coupling 
laser. Section IV investigates the role of the atomic mo- 
mentum. A conclusion completes our work in Section 
V. 

II. THREE-LEVEL SYSTEM AND LASER 
PROPAGATION 

A. Atoms 

We consider a cold sample of open three-level atoms 
interacting with two laser fields in the A configuration. 
The interaction scheme is shown in Fig. 0. The ground 
(or metastable) states |c) and |p), with energies respec- 
tively Ec and Ep, are excited to a common state |e), with 
energy E^, by two laser fields of frequencies ujc, i^p, elec- 
tric field amplitudes £c, £p, and wavevectors — i^a/c, 
with (a = c,p), c being light speed in vacuum. As stan- 
dard in laser without inversion and similar processes , 
the laser fields are indicated as coupling (c) and probe (p) 
lasers. The laser detunings are denoted by 6a — tOa—^ca 
and the Raman detuning from the two-photon resonance 
is Sr = Sc — Sp. The Rabi frequencies are given by 



with Vca the atomic dipole moment, supposed to be real 
for the sake of simplicity. The evolution of the atomic 
density matrix p is described by the OBE for an open 



system [|^. We assume the excited state relaxation, de- 
noted by A, due to spontaneous emission, as for a dilute 
atomic sample. In the open system the A rate is com- 
posed by the rates Fc and Fp for the decays into the 
ground states, and by the decay Fout into a sink state 
I out), not excited by the lasers 

A = Fe+Fp+Fout. (2) 

Thus the total decay rate of the excited state population, 
denoted by Fo, is equal to A, and the decay rate of the 
optical coherences, pcp and pcc, is equal to A/ 2. For a 
given open system the application of a repumping laser, 
as performed in , eliminates the decay channel out of 
the system. Thus we may compare the slow light prop- 
agation in an open system characterized by the above 
relaxation rates, to the propagation in the closed system 
produced by the repumping laser. That closed system is 
characterized by the following population decay for the 
excited state: 

rl = Te + Fp, (3) 

while in the presence of repumping the decay rate of the 
optical coherences remains A/ 2. In our analysis we will 
introduce also a decoherence rate 7cp for the ground state 
coherence pcp- 



le> 




lc> 

FIG. 1. Schematic representation of an open three-level 
atomic system interacting with two laser fields in the A con- 
figuration, with the decay rates Fc, Fp and Fout towards lower 
states. 

The response of the three-level system is determined 
by the presence of a noncoupled state, i.e. a state non- 
absorbing the laser radiation, given by 

|NC(i)) = i (l7p|c) -17ee^("---^)*b)) e-'^. (4) 
where fl is defined by 
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n = 
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The decoherence rate 7cp couples |NC(t)) to its orthog- 
onal state \C{t)). 

In both closed and open three-level systems an opti- 
cal pumping preparation of the noncoupled state, taking 
place within few excited state lifetimes, requires a signif- 
icant occupation of excited state population While 
in a closed system all the atomic population falls into 
the noncoupled state without loss of atomic population, 
for an open system a decay out of the three-level system 
takes place. For an open system, the occupation of the 
ground states |c) and \p) at long interaction times is dif- 
ferent from zero only at small Raman detunings, where 
coherent population trapping is efhcient. The analysis 
of Ref. 1 16 1 indicates that at long interaction times the 
slope of Re(pcp) at (5r reaches a constant nonzero 
value. In effect around (5r = the population is trapped 
in the nonabsorbing state up to an interaction time Q 
dependent on [|1^,0. A constant, and slow, group 
velocity is realized because the slope of the dispersion 
around 5r = does not change significantly for an in- 
creasing interaction time 0, even if the interval of sharp 
variation of the dispersion decreases with Q. However, 
the loss rate out of the three-level system increases with 
the Raman detuning, and the transparency window nar- 
rows for increasing interaction time. 



B. Lasers 

To describe the propagation of laser pulses through 
a cold sample of open three-level atoms, the electric 
field amplitudes, whence the Rabi frequencies, are as- 
sumed dependent on the time t and the coordinate z 
along the propagation direction. For slowly varying elec- 
tric field amplitudes, the Maxwell equations for the cou- 
pling/probe envelope Rabi frequencies expressed in terms 
of the pulse- localized coordinates z and t = t — z/c, 
rip(z,r) and f2c(z,r) reduce to jl^] 
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with a = (c, p) and where the parameter Kq is given by 

(7) 



N being the atomic density and eo the vacuum suscep- 
tibility. The Beer's absorption length for the probe field 
Zp in the absence of coupling laser is 



Cp 



A 
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In the ideal case, the solution of OBE and of Eq. (g) 
leads to a shape-invariant propagation of the probe pulse 



described through the dependence ^^{t — -^), Vg being 
the probe pulse group velocity. 

OBE for the density matrix and Eqs. (||) for the laser 
propagation were solved numerically for an initial Gaus- 
sian pulse on the probe transition 



npiz = o,t) = f^opcxp ( ^ 
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In the simulation the coupling laser, with constant Rabi 
frequency flc, was assumed switched on at earlier times 
t < 0, in order to realize a counterintuitive pulse se- 
quence. The atoms were initially prepared in the ground 
|p) state. Owing to the initial condition and the coun- 
terintuitive pulse sequence, the atoms were initially pre- 
pared in the noncoupled state. 

Because of the time dependence of the Rabi frequen- 
cies in the noncoupled state of Eq. (jj), the adiabatic 
coherent preparation requires a proper choice of the time- 
scales In STIRAP the conditions for the atoms 
in the initial \p) state to remain in the noncoupled state 
following the adiabatic evolution are 
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with a = {c,p). The first inequality, valid for resonant 
lasers, imposes a minimum Rabi frequency. Larger Rabi 
frequencies are required for detunings 5c and 5p different 
from zero [ p^ . The second inequality imposes a pulse 
width long compared to the inverse of the Rabi frequen- 
cies. We have verified numerically the role of both con- 
ditions on the slow light propagation. However it should 
be noted that the STIRAP conditions are required for 
an adiabatic transfer of the atomic population, while for 
slow light the atomic response is not important, and the 
only request is on the probe pulse propagation. For in- 
stance in the limit of a very small f7p Rabi frequency, a 
negligible STIRAP transfer is realized, still a slow light 
propagation may be produced. In Ref. [Q where the 
formstable pulse propagation was linked to the preserva- 
tion of population in the noncoupled state, the condition 
for a small nonadiabatic coupling between noncoupled 
|NC(t)) and coupled |C(<)) states was expressed through 
the following Rabi frequency coupling 



^ ^c^p 
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where the dot denotes the time derivative. On the basis 
of this Rabi coupling, the adiabaticity condition required 
for a formstable probe pulse propagation was written p^ ] 



(12) 



This relation, and the second one of Eq. (p^), should be 
satisfied for the slow light propagation. 

The probe Rabi frequency transmitted after the prop- 
agation through the cold gas medium was determined 
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numerically, with the pulse delay Td{z) at the given spa- 
tial position z defined by 

jT^r\n^iz,r)\'dr 



Tdiz) 



Finally the group velocity was calculated ml 
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(13) 
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The bandwidth Auj of the EIT transparency window 
determines the spectral components of the probe pulse 
propagating deep into the atomic medium without shape 
distortion. The EIT window was explicitly investigated 
in the applications of slow light to quantum entanglement 
and stopping of light |po| . The transparency window for 
both closed and open systems depends on the decay of 
the noncoupled state. For a closed system with = 
at (5r = the transparency window is given by 
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whence the window decreases towards a limiting value 
at large interaction times Q. For an open system that 
window is ilql 



Auj = 
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FIG. 2. Probe Rabi frequency ^p{z,t) as a function of 
the pulse-localized time r at diflFerent penetration distances 
in the conditions of ref. [Q, with Tc = A/3, Tp = A/2, 
and Tout = A/Q. The dashed line represents a reference 
probe pulse with no atoms in the propagation medium, scaled 
down by a factor two. The continuous and dashed-dotted 
lines are pulses for propagation through the medium at 
z — 30(p and at 2 = 63Cp. Continuous lines for jcp ~ 



and dashed-dotted lines for jcp = 10 



Initial Gaussian 



probe pulse width T = 80/ A, and initial coupling laser Rabi 
frequency fie ~ 0.18^1. 



For the probe pulse propagation, where the interaction 
time & can be approximated by the pulse duration T, in 
order to have a formstable pulse the transmission window 
should be larger than the Fourier frequency distribution 
of the pulse 



Auj > 1/T. 
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For an open system with bandwidth given by Eq. (p^), 
Eq. ( p^ ) imposes limitations on the Rabi frequencies. 
That condition is satisfied by the parameters of typical 
slow light pulses Q-^- Only for very short pulses, or 
very weak coupling lasers, the EIT bandwidth could pro- 
duce a large distortion of the pulse shape propagating 
through a three-level open system. We will derive in the 
next Section that, owing to Eqs. (^^ and (^2|), a bet- 
ter transmission of the slow light pulses is obtained by 
increasing the pulse duration T. 



III. PROBE TRANSMISSION 

Numerical results for the propagation of the probe 
pulse within a sample of open three level atoms are shown 
in Fig. ^. We used the parameters of the experiment in 
Ref. Q, i.e. Na D2 line parameters {A = 2tt ■ 5.9 MHz, 
A = 589.0 nm), cold atom density of 3.3-10^^ atoms/cm^, 
propagation length through the medium up to a distance 
z equal to 63 times the Beer's length Cp, initial Rabi fre- 
quencies flop and r^c around a tenth of A. The numerical 
results show a pulse peak height decreasing with the pen- 
etration distance, and a probe pulse propagating with a 
slow velocity. We calculated from the data of Fig. || the 
slow velocities shown in Fig. ||, i.e. in the range of the 
values measured in ref. Q. In order to compare these 
results quantitatively with the response of a closed sys- 
tem, the atomic parameters for the closed system should 
be chosen carefully. We may introduce a correspondence 
between closed and open systems by means of a repump- 
ing laser. As discussed previously, we compared the slow 
light propagation in an open system characterized by the 
Pee decay rate A of Eq. to that of a closed sys- 

tem characterized by the decay rate F^, of Eq. (||), the 
same optical coherence decay rate A/ 2 applying to both 
systems. Numerical analyses for a closed system corre- 
sponding to the open system of Fig. || have shown that 
the same pulse transmission and same slow light veloc- 
ity are obtained for the two systems. In fact the pulse 
propagation of Fig. ^ was examined for a counterintuitive 
coupling/probe pulse sequence and with laser parameters 
producing a very small excited state occupation, whence 
a small role of the Tout loss process. 

We have examined the probe transmission modifica- 
tion produced by the decoherence rate 7cp. The decoher- 
ence process decreases the |NC(i)) occupation, whence 
decreases the probe transmission. We have verified that 
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for 7cp larger than 0.001^, the probe transmission is well 
described through the EIT absorption length Cp ^""^ given 
by the following formula pi]]: 



EIT 



A 



27cpA 



(18) 
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FIG. 3. Probe-pulse velocity versus the inverse of the 
atomic density. The data are in (a) for coupling intensity 
of 12 mW/crn^ corresponding to coupling Rabi frequency Qc 
=0.56 A, and in (b) for a coupling intensity of 3 mW/cm^ 
corresponding to a coupling Rabi frequency Qc =0.18 A. The 
continuous line represents results for a simulation not includ- 
ing the atomic momentum, while results of the dashed line 
include in the simulation the modifications of atomic mo- 
mentum produced by the light forces. In (a) continuous and 
dashed lines cannot be distinguished on the scale of the plot. 

Note that Eq. (|l^) predicts no attenuation for the 
case of 7cp = 0. On the contrary the numerical results of 
Fig. ^ evidence an attenuation of the probe pulse even in 
the ideal case of zero ground state decoherence. In effect 
for the parameters of Fig. || the adiabaticity condition of 
Eq. ( |l^ required for the complete validity of the EIT de- 
scription is not fully satisfied. The adiabaticity condition 
was calculated for the laser pulses of Fig. H, as presented 
in Fig. U(a). strongly depends on the value of T, 
the probe pulse duration. Better adiabatic conditions, 
whence a larger probe pulse transmission, are realized at 
larger T values, with longer duration pulses, as shown in 
Fig. ^(b). We verified that the probe- light slow velocity 
has a weak dependence on the pulse length T, as in p6| . 
That weak dependence confirms that the dispersion slope 
of Re{pcp) versus the laser detuning reaches a constant 
value at large interaction time. 

We have examined the dependence of the slow light 
velocity on the atomic density n, in order to simulate the 
dependence on the cold atom temperature investigated in 



Q . In the experimental investigation the atomic density 
n was varied with the sample temperature. Our results 
for the slow light velocity are shown in Fig. ^ for two 
different values of the coupling laser parameters Oc. A 
linear dependence of the velocity on the inverse of the 
atomic density is derived from the data, with values be- 
tween few m/s and 100 m/s depending on the atomic 
density and the coupling laser Rabi frequency. The low- 
est velocities are reached by decreasing the coupling laser 
Rabi frequency fie. For the dense media required for 
cold atom slow light propagation, the reabsorption of the 
spontaneously emitted photons, neglected in our theoret- 
ical analysis, may play an important role. That reabsorp- 
tion is greatly reduced in the cigar shape geometries of 
cold atomic samples [p2|. 
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FIG. 4. In (a) non-adiabaticity test comparison between 
D,- and for T = 80 /A at 2 = 0. In (b) probe Rabi fre- 
quency Sip as a function of the pulse-retardation time t at 
z = 63("p in the open system of Fig. ^ with difi'erent tempo- 
ral widths of the probe pulse. Higher transmission pulses are 
obtained at T = 160/yl decreasing transmission at T = SO /A 
and r = 40/ A. 

The numerical simulation allowed us to verify the pres- 
ence of adiabatons on the coupling laser, i.e. modifica- 
tions in the coupling laser intensity propagating through 
the atomic sample synchronously with the probe laser 
psf . These modifications are small, less than a ten per- 
cent of the coupling laser intensity. Their observation 
in experiments as those of Refs. would represent a 

confirmation of the coupling/probe matched pulse prop- 
agation. 

The probe pulse propagation was examined also for the 
case of laser detunings i5c, i5p, different from zero, satisfy- 
ing the Raman resonance condition (5r = 0. The results 
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of such simulation are shown in Fig. |5|. For laser detun- 
ings different from zero, the probe field Rabi frequency 
acquires components in phase and out of phase with the 
initial probe laser field, and both components, Re{flp) 
and Im{flp) contribute to the pulse propagation. The 
absolute value |fip| is plotted in Fig. ||(a) while the sep- 
arate in-phase and out-phase components are shown in 
Fig. ^(b). The individual in-phase and out-phase com- 
ponents acquire complicate pulse shapes, even if the ab- 
solute Rabi frequency propagates preserving the initial 
shape, apart from an attenuation and a broadening. The 
presence of a laser detuning imposes stronger adiabaticity 
conditions, not well satisfied by the weak Rabi coupling 
frequencies required for slow light propagation. In fact, 
comparing the results of Fig. g(a) with those of Fig. ||, 
it appears that even for a small detuning the pulse shape 
of a slow pulse is less preserved. Increasing the laser 
detuning, the adiabaticity condition can be satisfied in- 
creasing the coupling laser Rabi frequency, at expenses of 
an increase of the slow velocity associated with the probe 
pulse. 




50 100 150 




50 100 150 

T ills) 

FIG. 5. Propagation of a probe pulse slow light for 
coupling/probe lasers detuned from the upper |e) level, 
maintaining the Raman resonance condition Sr = 0, i.e. 
Sc = Sp = —A, and other parameters as in Fig. ^. In (a) 
absolute value of the probe Rabi frequency, \flp\, versus the 
pulse-localized time r. As in Fig. ^the dashed line indicates 
a reference pulse. Continuous lines for different propagation 
distances through the medium. In (b) in-phase and out-phase 
components of the probe Rabi frequency, iie(Slp) and Im{Q,p) 
respectively, at propagation distance z = 63Cp. 

IV. ATOMIC MOMENTUM 

While the previous analysis neglected the atomic mo- 
tion, in a sample of cold atoms the modifications of the 
atomic momentum produced by the laser interactions 



should be taken into account. The atomic states will be 
classified by the internal variable \a) with (a = c,p,e) 
and by the external continuous variable p of the atomic 
momentum. For an atom initially in the \p) state with 
momentum p = 0, the Raman process of absorption 
and stimulated emission transfers the atoms from state 
\p,p — 0) into the state |c, fcp — fcc)- Including the atomic 
momentum the noncoupled state at t = becomes 

|NC,„(t = 0)) = ^ {n,\p, 0) - r!p|c, fcp - 4)) . (19) 

The kinetic energy Hamiltonian operator i/k = P^/2Af, 
M being the atomic mass, couples the |NCm) state to its 
orthogonal one |Cm) with a rate 7k given by 

= n ^ — 2M — (2^^ 

This rate 7k decreases the ground state atomic coher- 
ence, and whence modifies the slowing down process of 
the light pulse. For fee ~ fcp, a rate 7k equal zero is re- 
alized in the coupling/probe copropagating laser config- 
uration, and this result applies also for an initial atomic 
momentum p different from zero, as in the experiments 
of Refs. 1^,^. Instead the kc ^ —fcp counterpropagat- 
ing configuration was used in the laser cooling based on 
velocity-selective coherent population trapping For 
the experiment of |4| with orthogonal propagation direc- 
tions for coupling and probe lasers, and the parameters 
of Fig. H, the decoherence rate 7k is equal to 10^ s~^, 
close to the value used in the numerical analyses of Fig. 
^ with 7c ^ 0. It should be noted that the recoil fre- 
quency for sodium is wr, = hk'^/2M — 1.7 x 10^ s""'^, but 
in Eq. ( pO| ) the (ilp/fi)^ factor produces a smaller 7k. 

In order to test the influence of the atomic momentum 
decoherence on the slow light production, we solved nu- 
merically the generalized OBE, i.e. the density matrix 
equations including the atomic momentum p^ , for the 
parameters of the experiment in i.e. with orthogonal 
kc and fcp . The numerical results for the slow light veloc- 
ity including the atomic momentum in the propagation 
are shown in Fig. ^. The linear dependence of Vg on 1/n 
remains also when the atomic momentum is included in 
the analysis. The modification of the slow light velocity 
produced by the atomic momentum is around ten percent 
at the smaller coupling laser Rabi frequency. 

The modifications on the atomic momentum produced 
by the laser pulse interaction originate from the forces 
acting on the atoms. The longitudinal force, along the 
laser propagation direction, on the atom of the slow light 
medium contains both a radiation-pressure dissipative 
term and a dipole gradient reactive term ||2^. The dis- 
sipative radiation-pressure force F^p acts on the atoms 
for both resonant and non resonant probe lasers, while 
the reactive dipole force F^ip is present only for Sp ^ 0. 
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These forces are determined by the optical coherence of 
the atomic dipole matrix p6[: 



Frp{z, t) = -i-kpflpiz, t)pl^{z, t) + c.c. 
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where, owing to the counterintuitive pulse sequence, and 
the orthogonal kc , kp geometry of Ref . ||] , only the probe 
laser contributes to the force along the z axis. These 
forces, calculated for the propagating pulses of the pre- 
vious figures, are plotted versus time at z = in Fig. 
^ for the case of Sc — Sp = —A. From the figure it 
appears that for the chosen laser detuning the radiation 
pressure force is larger than the dipole force by few orders 
of magnitude. Similar values were obtained for the radi- 
ation pressure force for a probe laser in resonance, when 
the dipole force is equal zero. From the point of view of 
the atomic response, the total modification of the atomic 
momentum depends on the integral of the force over the 
interaction time, i.e. the pulse duration time. Because 
the radiation pressure force has an antisymmetric depen- 
dence on the time, its integral, whence the modification 
of the atomic momentum, is small. We have verified that 
for same laser parameters, detunings and Rabi frequen- 
cies, for instance those of Ref. |jl^, the time dependencies 
of the light forces are completely antisymmetric, so that 
the forces produce a total modification of the atomic mo- 
mentum equal zero. 




X(|1S) 

FIG. 6. Forces acting on sodium cold atoms versus the 
pulse-localized time t, at Sc = 5p — —A while the probe 
laser pulse of Fig. |^ enters the propagating medium (z=0). 
Radiation-pressure force frp in (a), and dipole force Fdip in 
(b). 



V. CONCLUSIONS 

We have examined the propagation of a probe pulse 
through a medium composed of an open three-level A 
system in the presence of a coupling laser acting on the 
adjacent transition of the system. Our numerical analy- 
sis confirms that also in an open system a large reduction 
of the probe-light velocity is obtained. The comparison 
with a closed system created by applying a repumping 
laser to the open system has demonstrated that for the 
parameters of the numerical analysis the velocities of the 
propagating pulse are equal in the compared open and 
closed systems. Also the amplitudes of the transmitted 
probe pulses are similar for the open and closed systems. 
This result corresponded to a set of laser parameters 
where the excited state occupation was very small, so 
that the external losses were not playing an important 
role on the atomic evolution. The intensity of the cou- 
pling and probe lasers were chosen in order to satisfy the 
requests of both EIT and excited state occupation. Fur- 
thermore a proper choice of the width of the probe pulse 
allowed us to satisfy the adiabaticity conditions of the 
atomic response, with a larger pulse transmission. We 
have examined the role of the laser pulse forces acting on 
the atoms. Those forces modify the pulse propagation 
and for the slowest pulses give an important contribution 
to the slow light velocity. The contribution of the atomic 
momentum to the slow light propagation can be tested 
using different laser propagation geometries, whence dif- 
ferent atomic momenta for the noncoupled state of the 
atomic preparation. 
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